Large Scale Magnetic Fields and the Number of Cosmic Ray Sources above 10^^ eV 
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We present numerical simulations for the two-point correlation function and the angular power 
spectrum of nucleons above 10^^ eV injected by a discrete distribution of sources following a simple 
approximation to the profile of the Local Supercluster. We develop a method to constrain the 
number of sources necessary to reproduce the observed sky distribution of ultra-high energy cosmic 
rays, as a function of the strength of the large scale cosmic magnetic fields in the Local Supercluster. 
While for fields B < 0.05/i G the Supercluster source distribution appears inconsistent with the data 
for any number of sources, fields of strength B ~ 0.3/i G could reproduce the observed data with a 
number of sources around 10. 



I. INTRODUCTION 

Despite a growing amount of data the origin of cosmic 
rays especially at the highest energies is still obscure. For 
ultra-high energy cosmic rays (UHECR) with energies 
above 10^^, there are still many open questions such as 
"How can particles be accelerated to these extremely high 
energies?" , and "What are their sources ?" ||l| . The best 
candidates for acceleration sources are powerful objects, 
such as hot spots of radio galaxies and active galactic 
nuclei but they are still not identified and it is still 
unknown how many of them contribute to the observed 
cosmic ray flux. 

The observed spectrum covers about 11 orders of mag- 
nitude, from 1 GeV to 10^^ GeV, and is described by 
a power law oc E~'^ with two breaks at the "knee" , at 
~ 4 X 10^^ eV, and at the "ankle", at 5 x 10^^ eV. Above 
the knee the spectrum steepens from a power law index 
7 ~ 2.7 to ~ 3.2. Above the ankle the spectrum flattens 
again to a power law index 7 ~ 2.8. Cosmic rays with en- 
ergies above the ankle cannot be confined by the Galactic 
magnetic field, and the lack of counterparts in our Galaxy 
suggests that the ankle marks a cross-over from a Galac- 
tic component to a component of extra-galactic origin. 
Data from the Fly's Eye experiment also suggest that 
the chemical composition is dominated by heavy nuclei 
up to the ankle and by protons beyond 

If UHECR have an extra-galactic origin, we would ex- 
pect a cutoff in the spectrum due to the fact that in 
the bottom-up scenario UHECR are assumed to be pro- 
tons accelerated in powerful astrophysical sources: Even 
if they can achieve, under extreme conditions, such high 
energies, they will lose their energy mostly by pion pro- 
duction on the microwave background. For sources fur- 
ther away than a few dozen Mpc this would predict a 
break in the cosmic ray flux known as Greisen-Zatsepin- 
Kuzmin (GZK) cutoff |], around 50EeV. This break has 
not been observed by experiments such as Fly's Eye 
Haverah Park §, Yakutsk §, Hires § and AGASA [|, 
which instead show an extension beyond the expected 
GZK cutoff and events above 100 EeV. 

One of the possible solutions to the lack of observed 
counterparts to the highest energy events P, 10 1 is to 



suppose the existence of large scale intervening magnetic 
fields with intensity B ^ 0.1 — 1 fiG ||l^, which would 
provide sufficient angular deflection even for high ener- 
gies and could explain the large scale isotropy of arrival 
directions observed by the AGASA experiment ||] as due 
to diffusion. 

It has been realized recently that magnetic fields 
as strong as ~ IfiG in sheets and filaments of large 
scale structures, such as our Local Supercluster, are 
compatible with existing upper limits on Faraday rota- 
tion 

In our previous paper fl^ we considered the effects 
of such strong magnetic fields in the particular case of 
a single source corresponding to Centaurus A, which is 
a radio-galaxy located in the southern hemisphere at a 
distance of 3.4 Mpc. There we employed detailed numer- 
ical simulations for the energy spectrum and the angu- 
lar distribution of ultra-high energy nucleons propagat- 
ing in extra-galactic magnetic fields of r.m.s. strength 
between 0.3 and 1 fiG. We found that this model is in- 
consistent with the data when B ~ 0.3/iG because the 
angular distribution predicted is not isotropic but con- 
centrated around the position of the source and because 
the northern hemisphere experiments should never have 
detected the highest energy events for which the angular 
deflection is too weak to bring the particle in the field of 
view of these experiments; therefore we argued that at 
least a few sources within the GZK cutoff are required to 
produce the observed UHECR flux. 

The goal of our present paper is to elaborate more 
detailed constraints on the number of sources necessary 
to reproduce the observed distribution, as a function of 
the poorly known strength of the extra-galactic magnetic 
field in our Local Supercluster. As will be explained be- 
low in more detail, we assume a discrete distribution of 
sources in the Local Supercluster permeated by magnetic 
fields of strength up to i? ~ 0.3/i G. 

As in our previous paper, we restrict ourselves to 
UHECR nucleons, and we neglect the Galactic contribu- 
tion to the deflection of UHECR nucleons since typical 
proton deflection angles in galactic magnetic fields of sev- 
eral fiG are < 10° above 4 x lO^^ eV [|l|, ^ , and thus are 
small compared to deflection in > 0.3/^G fields extended 
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over megaparsec scales. 



II. NUMERICAL SIMULATIONS 



As statistical quantities used to test various scenarios 
we adopt the angular power spectrum based on the set of 
spherical harmonics coefficients aim , as used in Ref . [ pT[ , 
which is sensitive to anisotropics on large scales, and 
the two-point correlation function as defined in Ref. , 
which contains information on the small scale anisotropy. 

As will become apparent, the statistics for these quan- 
tities is so far limited by the small number of observed 
events but the present development of large new detec- 
tors will considerably decrease their statistical uncertain- 
ties. In particular, the Pierre Auger experiment [ p^ will 
combine ground arrays measuring lateral shower cross 
sections with fluorescence telescopes measuring the lon- 
gitudinal shower development. Since two of these hy- 
brid detectors are planned, one in the southern hemi- 
sphere currently under construction in Argentina, and 
one in the northern hemisphere, full sky coverage will be 
achieved, with an exposure that is practically uniform in 
right ascension, and a geometrical dependence on decli- 
nation. There are furthermore plans for space based air 
shower detectors such as OWL a-^d EUSO |21 which 
may also achieve full sky coverage. For this reason it is 
feasable to consider a multi-pole analysis of the angular 
distribution which involves statistical estimators of inte- 
grals covering the full sky: Since these estimators involve 
factors l/iOi, where uji is the exposure associated with 
the ith observed direction, they are undefined if the ex- 
posure vanishes anywhere on the sky. However, even in 
the absence of full-sky coverage one can define analogous 
quantities and their estimators (which are then different 
from the usual spherical multipoles) by simply restricting 
them to the area of the sky where the exposure function 
does not vanish. We will use these modified statistical 
quantities to compare model predictions with the exist- 
ing AGASA data. 

The auto-correlation analysis provides information 
about the small scale anisotropy and can be applied to 
partial sky coverage such as for the AGASA experiment 
without restriction. The observed data actually show sig- 
nificant small-scale angular clustering (five doublets and 
one triplet within 2.5° out of 57 events above 40 EeV). 
This clustering has a chance probability of less than 1% in 
the case of an isotropic distribution. It has been pointed 
out that in the presence of turbulent extra-galactic mag- 
netic fields of fractions of a micro Gauss clustering could 
be induced by magnetic lensing |2|, |2^, The 
auto-correlation analysis presented in the present paper 
will demonstrate that quantitatively. 

The paper is organized as follows: in section II we 
briefly describe our numerical simulations, in sections III 
and IV we present our results on multi-pole analysis and 
auto-correlation function, respectively. Section V briefly 
reconsiders Centaurus A as the unique source in case of 
a field as strong as a micro Gauss and in section VI we 
conclude. 



We use the same numerical approach used in earlier 
publications |2j], but we take a discrete distri- 

bution of sources centered at 20 Mpc from Earth and 
distributed on a sheet with a Gaussian profile of thick- 
ness 3 Mpc and radius 20 Mpc, with both magnetic field 
strength and source density following the profile of the 
sheet and no sources present within 2 Mpc from the ob- 
server. This is a simple approximation to our location in 
the Local Supercluster and to its shape. We also assume 
that the sources inject protons with a E~^ '^ spectrum 
extending up to ~ 10^^ eV. We note that the angular dis- 
tributions are not very sensitive to assumptions on the 
injection spectrum. 

We assume a random turbulent magnetic field with 
power spectrum (i3(fc)^) cx /c"^ for 2tt / L < k < 2tt/Ic 
and {B'^ik)) = otherwise. We use ub = —11/3, corre- 
sponding to Kolmogorov turbulence, in which case L, the 
largest eddy size, characterizes the coherence length of 
the magnetic field. For the latter we use L ~ 1 Mpc, cor- 
responding to about one turn-around in a Hubble time. 
Physically one expects Ic ^ L, but numerical resolution 
limits us to 1^ > 0.008L. We use Ic ^ 0.01 Mpc. The 
magnetic field modes are dialed on a grid in momentum 
space according to this spectrum with random phases and 
then Fourier transformed onto the corresponding grid 
in location space. The r.m.s. strength B is given by 
S2 = /g" dkk'' {B\k)). 

Typically, 5000 trajectories are computed for each re- 
alization of the magnetic field obtained in this way and 
of the source positions, for 10-20 realizations in total. 
Only those trajectories that cross a sphere of 1.75 Mpc 
radius around Earth (corresponding to 5° viewed from 
20 Mpc distance) are used. Each time such a trajectory 
crosses this sphere, arrival direction and energy are reg- 
istered as one event. Each trajectory is followed for a 
maximal time of 10 Gyr and as long as the distance from 
the observer is smaller than double the source distance. 
The results do not significantly depend on these cut-offs. 
Furthermore, the distance limit is reasonable physically 
as it mimics a magnetic field concentrated in the large 
scale structure, with much smaller values in the voids, as 
generally expected. Similar codes have been developped 
in Rcfs. (2§1. 

When dialing simulated data sets from the simulated 
sky distributions, one has to take into account the non- 
uniform exposure of the particular experiment consid- 
ered. This can be done by dialing from the simulated 
distribution multiplied by an exposure function depend- 
ing on the sky solid angle fi. This function, measured 
in units of km^years, gives the effective time-integrated 
collective area of the detector in a given direction Q.. A 
detector which operates continuously will have an expo- 
sure function roughly independent of right ascension and 
thus will only depend on the declination angle 5. We will 
only need the exposure function up to an irrelevant over- 
all normalization. For a detector at a single site we use 



3 



the following parameterization: 

w{S) oc cos oq cos S sin am + ctm sin oq sin S , 



(1) 



where oq is the latitude of the detector and am is zero for 
^ > 1, TT for ^ < —1, and cos~^(^) otherwise, where ^ = 
(cos^m — sin Co sin (5)/ [cos oq cos 5]. The angle dm is the 
maximal zenith angle out to which the detector is fully 
efficient (60° for Auger, 45° for AGASA). The exposure 
function which results for the AGASA experiment and 
that we will use in the following, has been discussed, for 
example, in Ref. see in particular Fig. 2 there. 



III. THE ANGULAR POWER SPECTRUM 



The angular power spectrum is defined as the average 



C{1) 
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(2) 



and the statistical estimator for the spherical harmonic 
coefficients aim is given by |]l7| 



1^1 

azm = T7 —Yimiu") , 
A/ ^ 

i— 1 



(3) 



where N is the number of discrete arrival directions, ei- 
ther of the real data, or randomly sampled from the simu- 
lated sky distributions. Furthermore, u>i the total exper- 
imental exposure at arrival direction u^, M — X^iLi 1/^* 
the sum of the weights and Yim is the real- valued 

spherical harmonics function. 

In order to obtain the statistical distribution of the 
C{1) predicted by specific simulated scenarios, we dial 
C{1) typically lO"' times from the simulated distributions 
[multiplied by the exposure function uj{5)\ for each real- 
ization of the magnetic field and the source positions. 

For each I we plot the average over all trials and re- 
alizations as well as two error bars. The smaller error 
bar (shown to the left of the average) is the statistical 
error, i.e. the fluctuations due to the finite number N 
of observed events, averaged over all realizations, while 
the larger error bar (shown to the right of the average) is 
the "total error", i.e. the statistical error plus the cosmic 
variance, in other words, the fluctuations due to finite 
number of events and the variation between different re- 
alizations of the magnetic field and source positions. 

To estimate the true power spectrum from Eqs. (|^) 
and (|^) requires data with full sky coverage and therefore 
at least two detector sites such as forseen for the Auger 
experiment. For its exposure function wc add Eq. for 
two sites located at ao = —35° and at ao ~ 39°. The 
AGASA experiment only has partial sky coverage and, 
consequently, the true multi-pole spectrum cannot be 
computed from its data. For this case we consider the 
quantities defined by restricting Eq. (H) to the sky area 



where uj{5) > 0. This method is also used in the analysis 
of cosmic microwave background fluctuations where win- 
dow functions are used which are unity in the observed 
region and zero elsewhere. In our case this corresponds to 
using the AGASA exposure function for the uoi in Eq. (|^) 
for the coefficients aim- This defines the modified angu- 
lar power spectrum C{1) both for the simulated data sets 
and the real data. 

We start by comparing in Fig. ^ the power spectra 
predicted by the completely isotropic distribution with 
the AGASA exposure function with the actual AGASA 
results which appear completely consistent with isotropy 
on large scales. Note that the increasing power for I = 
and / = 1 is due to the incomplete sky coverage of 
AGASA. 



58 events with 
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FIG. 1: Comparison of the angular power spectrum C(Z), 
Eqs. (P), (||), resulting from the A'' = 58 events above 40 EeV 
observed by AGASA (histogram), with the one predicted for 
an isotropic distribution (diamonds with error bars represent- 
ing the statistical error), as a function of multipole I. 



In the following figures, in case of full sky coverage, 
we show as solid line the analytical prediction for an 
isotropic distribution. In this case the power is the same 
for all Z-values and decreases as l/N as the number of 
arrival direction increases. AGASA data are shown as 
histograms. A pure mono-pole intensity distribution is 
equivalent to isotropy while the strength of other multi- 
poles relative to the mono-pole is a measure of anisotropy. 

Since the typical experimental angular resolution is ~ 
3°, in principle information is contained in modes up to 
I ~ 60. In the following we show the values of C{1) with 
I only up to 10 because the structure on small scales 
corresponding to larger / is better described by the auto- 
correlation function discussed in the next section. 

In Figs. ^ and H, we compare the angular power spec- 
trum C{1) predicted for the AGASA experiment at en- 
ergies E > AOEeV, for magnetic field strength B = 
0.05/xG, and 100 and 400 sources, respectively, with the 
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actual AGASA data. Both plots have been obtained for 
iV = 58, the present number of events with energies above 
4 X 1019 eV observed by AGASA. 
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percluster: For a number of sources much larger than 
the number of observed events N, it is likely that each 
observed event has been produced by a different source. 
The number of contributing sources is thus maximal and 
the fluctuations around the assumed (non-isotropic) dis- 
tribution is minimal, making the anisotropy more visible. 
This is illustrated by Fig. H, which shows the UHECR 
angular distribution as seen on Earth in terrestrial coor- 
dinates for E > 40EeV, B = O.OSAiG, and 400 sources. 
The distribution is concentrated around the solid line 
which represents the Supergalactic plane. 
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FIG. 2: The angular power spectrum C(l), Eqs. (|2|), (|3|), as 
a function of multipole I, obtained for the AGASA exposure 
function, for A*' = 58 events observed above 40 EeV, sam- 
pled from 10 simulated realizations for B = 0.05/iG with 100 
sources in the Local Supercluster. The diamonds indicate the 
realization average, and the left and right error bars repre- 
sent the statistical and total (including cosmic variance due 
to different realizations) error, respectively, see text for expla- 
nations. The histogram represents the AGASA data. 
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FIG. 3: Same as Fig. ^, but for 400 sources. 

The case of 100 sources seems roughly consistent with 
the experimental data while the case of 400 sources shows 
some deviations for the lowest multi-poles. This can be 
interpreted as the magnetic field being too weak to suf- 
ficiently isotropize the arrival directions with respect to 
the sources which were assumed to follow the Local Su- 



FIG. 4: The angular image in terrestrial coordinates, averaged 
over all 10 magnetic field realizations of 5000 trajectories each, 
for events above 40 EeV, as seen by a detector covering all 
Earth with B — 0.05fi G and 400 sources. The grey scale 
represents the integral flux per solid angle. The solid line 
marks the supergalactic plane. The pixel size is 1° and the 
image has been convolved to an angular resolution of 2.4° 
corresponding to AGASA. 

In contrast, the scenario with 100 sources seems to be 
more sensitive to the limited statistics due to the rela- 
tively small number of events observed by AGASA. In 
fact, the statistical errors due to the small number of 
events at low multi-poles is higher in Fig. || than in Fig. ||. 

For full sky exposure function corresponding to the 
Auger parameters, and assuming 500 oberved events, we 
obtain the situation shown in Fig. ||. In this case the 
deviation from isotropy, plotted as the solid line, is much 
more evident and would be easily determined by future 
observations. As we will show and explain in the next 
section, the scenario with 100 sources is, however, ruled 
out by the AGASA data from the analysis of the auto- 
correlation function. For the relatively small deflection 
induced by B = 0.05^ G, the number of sources must 
be at least as large as the number of events observed 
in different directions; much fewer than 100 sources are 
therefore ruled out in this case. 

We now investigate whether stronger magnetic fields, 
by providing larger angular deflection, might provide a 
better match to isotropy. In particular, we focus on the 
case where B = 0.3^ G. In Fig. | we show results for 
B = 0.3 /iG and 10 sources, all other assumptions being 
the same as in Fig. |^. 

The distribution seems to be roughly consistent with 
the data, but we also found that 5 and 100 sources result 
in almost the same distribution. Since in this case the 
limited statistics does not allow us to discriminate be- 
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FIG. 5: The angular power spectrum C{1), Eqs. (Q), as 
a function of multipole I, obtained for the Auger exposure 
function, assuming TV = 500 events observed above 40 EeV, 
sampled from 10 simulated realizations for B = 0.05/iG with 
100 sources in the Local Supercluster. Average and error 
bars are as in Fig. ^. The solid line represents the analytical 
prediction for an isotropic distribution. 
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FIG. 6: Same as Fig. ^, but for B = 0.3^ G and 10 sources. 



tween widely different number of sources, we turn to the 
case of Auger exposure with full sky coverage, assuming 
500 events observed above 40 EeV. The results are shown 
in Fig. and Fig. || for the case for 5 and 10 sources, re- 
spectively. The case of 100 sources is already ruled by 
the auto-correlation function of the AGASA data, as will 
be shown in the next section. 

Note that as for the weak field case the scenarios shown 
in Fig. and ^ predict an anisotropy that should be easily 
detectable by the Pierre Auger experiment, in contrast 
to the AGASA experiment (compare Fig. ||). More gen- 
erally, we note that the scenario with B = 0.05/iG gives a 



FIG. 7: Same as Fig. ^ but for B — 0.3/i G, 5 sources, and 
19 realizations. 
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FIG. 8: Same as Fig. H, but for B = O.Sfi G and 10 sources. 



distribution very different from the one with B — O.SfiG. 
Thus future experiments with full sky coverage should 
be able to give important information on the strength of 
the magnetic field in the Local Supercluster by perform- 
ing a multi-pole analysis. Furthermore, in the present 
case of i? = O.Sfi G Figs. ^ and || show that multipoles 
I ^ 3 hardly depend on the number of sources, whereas 
the lowest multipoles , which are less influenced by de- 
flection, have a noticeable dependence on the number of 
sources, similarly to the case of weak field B < 0.05/x G. 

If a nearly isotropic angular distribution is confirmed 
by future observations, we can conclude that for magnetic 
fields B ~ 0.3/iG a number of sources of 5-10 would 
be favored by the data. As will be shown in the next 
section, the auto-correlation function does not allow a 
much higher number of sources because magnetic lensing 
would not produce sufficient clustering on small scales. 
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On the other hand, much fewer than 5 sources are ruled 
out by the arguments given in our previous paper . 

Assuming an energy independent exposure function 
and using a simple spectrum, it is possible to es- 
timate the number of events which will be observed in 
the future. In the case of Auger observatories, for a to- 
tal acceptance of ~ 7000 km^ sr per array, in five years 
we should observe ^ 2200 events above 4 x 10^^ eV 
Here we have been conservative and used N — 500. 



IV. AUTO-CORRELATION FUNCTION 

For the auto-correlation function we follow the same 
approach used in Ref. |]l8j. We start from either actual 
data or from a randomly generated set of N events di- 
aled from the simulated distributions, multiplied by the 
exposure function. For each event we divide the sphere 
into concentric bins with a fixed angular size A9, and 
we count the number of events falling into each bin. We 
then divide by twice the solid angle size S{9) of the cor- 
responding bin, arriving at 



N{0) 



1 



2S{0) 



Ri] (6') : 



(4) 



where 



Rij (6) 



1 if 9ij is in same bin as 9 
otherwise 



We note that the auto-correlation function in the strict 
sense would include a factor (A/'^WiWj) under the sum 
in Eq. However, the differences are small and in 

any case we are free to choose any statistical quantity as 
long as the same quantity and its fiuctuations are used 
to compare simulations and data. 

In analogy to the previous section, for each magnetic 
field and source position realization we dial N{9) 
times from the simulated distributions in order to obtain 
its average and variances for which wc plot the same two 
error bars as for the power spectrum. The histograms 
shown subsequently represent again the result for the 
AGASA data, where the sharp peak at small separation 
angles results from the six observed clusters. We have 
verified that using incorrect exposure functions in gen- 
eral destroys the agreement found between a simulated 
isotropic distribution and the data at large 9 [|j . The ob- 
served distribution of events thus reflects the non-uniform 
exposure H, |2|. 

We start by comparing the autocorrelation function for 
the real AGASA data with the isotropic distribution in 
Fig. ||. This demonstrates that the AGASA data are com- 
pletely consistent with isotropy except at scales larger 
than a few degrees. 

In Fig. |lo|we show the angular correlation function for 
iV = 58 events with energies E > AOEeV, predicted by 
simulations with B = 0.05/^G and 100 sources, using the 
AGASA exposure function. 
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FIG. 9: Comparison of the angular correlation function N{6), 
Eq. (^), resulting from the N = 58 events above 40 EeV ob- 
served by AGASA (histogram), with the one predicted for an 
isotropic distribution (diamonds with error bars representing 
the statistical error), as a function of angular distance 0. A 
bin size AO = 2° was used. 



58 events with 40. < E/EeV < 1000. 



12 



10 



8 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



I I I I I I I I I I 




I I 



I I I I I 



10 20 30 40 

theta (deg) 



50 



60 



FIG. 10: The angular correlation function N{6), Eq. (^, as 
a function of angular distance 6, obtained for the AGASA 
exposure function, for A*' = 58 events observed above 40 EeV, 
sampled from 10 simulated realizations for B = 0.05/.iG with 
100 sources in the Local Supercluster. Average and error bars 
are as in Fig. ^. The histogram again represents the AGASA 
data. A bin size A9 = 2° was used. 



This case shows no correlation at angles as small as the 
angular resolution, where AGASA shows a peak, whereas 
there are strong correlations at larger angles, which is 
not consistent with the observed isotropic distribution at 
large scales. This also corresponds to the fact that in the 
case of weak magnetic fields we expect that clusters just 
refiect the point-like sources but if the number of sources 
is much larger than the number of observed events N, 
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each source constributes at most one events and cluster- 
ing is not possible. As remarked in the previous section, 
a much smaller number of sources is not possible either 
due to the large number of observed arrival directions. 

In Figs. |l^, and ^we show the angular correlation 
fimctions predicted by scenarios with B = 0.3/iG, with 
5 and 10 and 100 sources, respectively. In the case of 
100 sources the simulated distribution do not show any 
correlation at small scales. Similarly to the weak field 
case, this can be understood due to the fact that the 
source images produced by magnetic lensing contain at 
most one event if the number if sources is much larger 
than N. Note also that cosmic variance becomes very 
small for N > 100. 



58 events with 40. < E/EeV < 1000, 



10 

_ 8 
'a 
o 

5 6 

z 

-1 



I I I I I I I I I I I I I I 





I I 



I I I I I 



10 20 30 40 

theta (dcg) 



50 60 



FIG. 11: Same as Fig. but for B = O.SfiG and 5 sources 
with 20 realizations. 
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FIG. 13: Same as Fig. |l|, but for B = 0.3^i G and 100 sources 
with 18 realizations. 



imate current upper limit for the number of sources. 
On the other hand for 5-10 sources the simulated auto- 
correlation function seems to be in agreement with the 
observed clustering at small scales. In Fig. |l^ we show for 
one of these cases what could be expected for Auger ex- 
posures. This demonstrates that, for the amount of data 
expected with next generation experiments, the statis- 
tics will be dominated by cosmic variance instead of the 
limited number of events observed, as is presently the 
case. 
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FIG. 12: Same as Fig. 0, but for B = 0.3^ G and 10 sources. FIG. 14: Same as Fig. |l|, but obtained for the Auger exposure 

function, assuming TV — 500 events observed above 40 EeV. 

We obtain the same result for B = O.l/xG and 100 
sources. Thus we can argue that 100 is an approx- 
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V. ONLY ONE SOURCE: CENTAURUS A 

Now we brieflyreconsider the model discussed in our 
previous paper [|4|, with Centaurus A as single source, 
and B = IfiG. The predictions for the angular power 
spectrum and the auto-correlation function are compared 
with the AGASA data in Figs. |l^ and |l^, respectively. 
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FIG. 15: Same as Fig. ^, but for the scenario with Centaurus 
A at 3.4 Mpc distance as the single source, a field of 1/iG, 
and 20 realizations. 



images the source are produced ||2^. Correlations up 
to relatively large scales would only appear at energies 
above ~ 10^° eV, as can be seen in Fig. |l^ which was 
produced for the Auger exposure function. 



events with ICQ. < E/EeV < 1000. 
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FIG. 17: The angular correlation function N{6) as a func- 
tion of angular distance 6, obtained for the Auger exposure 
function, assuming A'^ = 68 events observed above 100 EeV, 
for the scenario with Centaurus A at 3.4 Mpc distance as the 
single source, a field of 1/i G, and 20 realizations. The bin size 
is again AS = 2°. 
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FIG. 16: Same as Fig. but for the scenario with Centaurus 
A at 3.4 Mpc distance as the single source, a field of l/^G, 
and 20 realizations. 

The angular power spectrum shows a ~ 3cr deviation 
from the data dXl — 2. Furthermore, the auto-correlation 
function does not show significant correlations at angular 
resolution scales, in contradiction to the data. This is due 
to the fact that for a magnetic field as strong as 1/iG, 
we are in a range of energies where many overlapping 



VI. DISCUSSION AND CONCLUSIONS 

In the present work we assumed a discrete distribution 
of sources distributed in the Local Supercluster and esti- 
mated the number of sources necessary to reproduce the 
experimental data, in dependence on the typical strength 
of the extra-galactic magnetic fields permeating the Local 
Supercluster. As statistical quantities for this analysis we 
used spherical multi-poles and the auto-correlation func- 
tion. We found that for weak magnetic fields ^ 0.05/iG 
the simulation predictions appear to be not consistent 
with the observed distribution for any number of sources 
because the magnetic field is too weak to isotropize the 
anisotropic distribution associated with the Supcrgalac- 
tic plane. Full sky experiments of the size of the Pierre 
Auger project will be sensitive to the difference in the dis- 
tribution of multi-poles between weak and strong mag- 
netic fields, which thereby could give direct informations 
about the strength of the magnetic fields. For stronger 
magnetic fields ~ 0.3/iG we found that the number of 
sources is constrained. In our previous paper JT^ we al- 
ready showed that a single source cannot reproduce the 
observed isotropic distribution. Here we found that for 
> 100 sources the auto-correlation function does not re- 
produce the correlations observed at small scales. This 
can be interpreted by the fact that for a number of 
sources much higher than the number of observed events 
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there are more source images produced by lensing than 
observed events and thus clustering is not observable. 
Therefore, the current upper limit on the number of con- 
tributing sources is ~ 100 in this case. On the other 
hand, a number of sources around 10 seems to repro- 
duce quite well the observed small scale clustering. We 
also showed that the model with Centaurus A as the 
only source and very strong fields ~ l/^G, considered 
as marginally consistent in our previous paper |Q, is 
not consistent with isotropy at large scales due to a pre- 
dicted quadrupole deviation from isotropy and because 
the auto-correlation function is not consistent with the 
clustering at small scales observed by AG AS A. We con- 
clude then that a distribution of ~ 10 sources in the Local 
Supercluster, with magnetic fields in the sub micro Gauss 
range, could reproduce at least current observations. Our 
approach can equally be applied to other source distri- 
butions. 

Due to the still sparse statistics of current data, in 
the present paper we refrained from quantifying the sta- 
tistical significance of deviations between models and 
data, because small number fluctuations are in general 
not Gaussian and, for different multipoles and separation 



angles, can be correlated. Results based on comparison 
with the AGASA data presented here should rather be 
understood as suggestive tendencies. Quantitative signif- 
icances can be obtained by determining in how many sim- 
ulated trials a certain quantity, such as the multi-poles 
and auto-correlations studied here, or certain combina- 
tions thereof, show deviations from the data of opposite 
sign to the average deviation. We leave that to a study 
with data of much higher statistics above ~ 10^^ eV than 
available today, as expected from future full-sky experi- 
ments. These experiments will put much more stringent 
constraints both on the number of sources and the mag- 
netic field strength. 
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